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ABSTRACT 

Context. Plasma processes close to supernova remnant shocks result in the amplification of magnetic fields and in the acceleration of 
electrons, injecting them into the diffusive acceleration mechanism. 

Aims. The acceleration of electrons and the magnetic field amplification by the collision of two plasma clouds, each consisting of 
electrons and ions, at a speed of 0.5c is investigated. A quasi-parallel guiding magnetic field, a cloud density ratio of 10 and a plasma 
temperature of 25 keV are considered. 

Methods. A relativistic and electromagnetic particle-in-cell simulation models the plasma in two spatial dimensions employing an 
ion-to-electron mass ratio of 400. 

Results. A quasi-planar shock forms at the front of the dense plasma cloud. It is mediated by a circularly left-hand polarized elec- 
tromagnetic wave with an electric field component along the guiding magnetic field. Its propagation direction is close to that of the 
guiding field and orthogonal to the collision boundary. It has a frequency too low to be determined during the simulation time and a 
wavelength that equals several times the ion inertial length. These properties would be indicative of a dispersive Alfven wave close to 
the ion cyclotron resonance frequency of the left-handed mode, known as the ion whistler, provided that the frequency is appropriate. 
However, it moves with the super-alfvenic plasma collision speed, suggesting that it is an Alfven precursor or a nonlinear MHD wave 
such as a Short Large-Amplitude Magnetic Structure (SLAMS). The growth of the magnetic amplitude of this wave to values well in 
excess of those of the quasi-parallel guiding field and of the filamentation modes results in a quasi-perpendicular shock. We present 
evidence for the instability of this mode to a four wave interaction. The waves developing upstream of the dense cloud give rise to 
electron acceleration ahead of the collision boundary. Energy equipartition between the ions and the electrons is established at the 
shock and the electrons are accelerated to relativistic speeds. 

Conclusions. The magnetic fields in the foreshock of supernova remnant shocks can be amplified substantially and electrons can be 
injected into the diffusive acceleration, if strongly magnetised plasma subshells are present in the foreshock, with velocities an order 
of magnitude faster than the main shell. 

Key words, acceleration of particles, methods: - numerical, plasmas 

1. Introduction the shock transition layer repeatedly. Electrostatic instabilities 

dominate for nonrela tivistic flows in unmagnetized plasmas 

Supernova remnants (SNRs) emanate energetic electromagnetic jBret et ani2008t lBretll2009i) and they can neither accelerate the 

radiation, which demon strates the acceleration of electrons to electrons to highly relativistic speeds dSircombe et al.ll2006l) nor 

ultrar elativistic speeds ('Tanim ori et al.l 119981: lUchivama et al.l amplify the magnetic fields. 

120071) and t he generation or amphfication o f magnetic fields j., „i„„frnn>i muM he arreleratpH hv nla^sma based 

X. \noHr^;.^,Jonnell\7„iuTf oi lhnn«:K Th» electrons could be acc elerated by p lasma based 



^Elhson & yiadimrro 3200i| Volk et aL lM). The Itkely ongm ^j^^ ^^.^j^ accelerators teinghametall by 

f 'cxTo 'f l/ , ° onn^ Strong magnetic fiel ds is electron surfing accelera tion (iKatsouleas & DawsonI 

the SNR shock (lMarcowiflietal. ,.2006; PelletieiiliDIIOOi. bhira & Takaharal l2007t IPiecLann et alJ l2008al). double 



T he nonrelativistic expansion speed of the m ain SNR blast layers jRaadu & RasmussenI Il988t iDieckmann & Bre3 l2009l) 

shell dKulkarni et al.l [19981 iFransson et all '2002*) and the weak or by processes that exploit a velocity shear in the plasma 

magnetic field of the ambient medium ([Volk et al. 2005), into outflow (Rieger & Duffy 2006). It has, however, not yet been 

which this shell is expanding, are obstacles to the magnetic field demonstrated with multi-dimensional and self-consistent simu- 

amplification by plasma instabilities and to the electron accel- lations that these mechanisms can indeed achieve the required 

eration out of the thermal plasma distribution to moderately electron acceleration and magnetic field amplification. The 

relativistic energies. Such an acceleration is needed for their non-relativistic shocks, which are found between the main SNR 

injection (Car gill & Papadopoulos 1988; Kirk & Dendy 2001; blast shell and the ambient medium, can also probably not 

iKuramitsu & Krasnosel skikh 2005a b) into the diffusive shock transfer significant energy from the ions to the electrons and 

acceleration process (See iDrury. (,1983i) ) so that they can cross accelerate the latter to relativistic speeds dHoshino & Shimadal 
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2004t IScholer & Matsukivol l2004t 



Lembege et alj 

Chapman et al.' 2005; 'Sor asio et al.l 120061 : lAmano & Hoshinol 



2007; Umeda et al. 2009). 

A viable acceleration mechanism may develop ahead of 
the main SNR shock, if we find subshells that outrun the 
main blast shell. These subshells may move faster than the 
typical peak speed of 0.2c of the main shell. An expansion 
speed as high as 0.9c may have been observed for a sub - 
shell ejected by the supernova SNl 998b w ( Kulkarni et aljl998l) . 
Most supernovae are less violent and their subshells are prob- 
ably slower The density of the subshell plasma is well be- 
low that of the main blast shell and its dynamics will be in- 
fluenced to a larger extent by the upstream magnetic field than 
the dynamics of the latter This is true in particular, if the 
upstream magnetic field has been pre-amplified by the cos- 
mic rays (Winske & Leroy 1984; Bell 2004; Reville et alJl2007t 
iNiemiec et al.l2008l:IRiquelme & Spitkovskyl2009l) . A fast mag- 
netized shock would form in the foreshock of the main SNR 
shock that can result in a stronger magnetic field amplification 
and electron acceleration. 

We examine with a particle-in-cell (PIC) simulation the for- 
mation of a shock in a plasma, in which a strong guiding mag- 
netic field is quasi-parallel to the plasma expansion direction. 
Whistler waves, which become Alfven waves at low frequen- 
cies, occuring at such sho cks can be effi cient electron accel- 
erators (Levi nson 1992; Kuramitsu & Kra snoselskikh 2005a b; 
iMiteva & Man n 2007; McClements & Fletc her 2009). Whisders 
and Alfven waves are circulary polarized if they propagate par- 
allel to the guiding magnetic field. We briefly summarize their 
properties and focus on the low-frequency modes with a left- 
hand polarization. These modes are qualitatively similar to the 
quasi-parallel propagating ones we consider here. A more thor- 
ough description of the dispersion relation of two-fluid waves 
and the shift of the resonance fre quencies for quasi-parallel 
prop agation can be found elsewhere (ITreumann & BaumiohannI 
[1997.) . 

The dominant waves, which we will observe, grow in a 
plasma with an electron gyrofrequency that exceeds the plasma 
frequency. As we increase the frequency in such a strongly mag- 
netized plasma towards the ion cyclotron frequency, the Alfven 
waves with a left-hand circular polarization become dispersive. 
These waves resonate with the ions and their frequencies remain 
below the ion cyclotron frequency. Alfven modes with a left- 
hand circular polarization just below the ion cyclotron frequency 
are called ion whistlers. Whistlers are predominantly electro- 
magnetic for small propagation angles relative to the guiding 
ma gnetic field, as it has been discussed for high-frequency ones 
by jTokar & Garvlll985h . and if they have low wavenumbers. 
Ion whistlers or dispersive Alfven waves develop a field-aligned 
electric field component close to the resonance frequency, by 
which they can inte ract nonhnearly with the plasma particles 
and accelerate them jMcClements & Fletcher 2009). Any wave 
growth will furthermore result in an increasing energy density of 
the magnetic field. 

Alfven waves and other magnetohydrodynamic (MHD) 
waves can grow to amplitudes, at which they start to interact 
nonlinearly with the plasma (Stasiewicz & Ekeberg 2008). Short 
Large Amplitude Magnetic Field Structures (SLAMS) are non- 
linear MHD waves occuring at quasi-parallel shocks and may be 
relevant for our si mulation. The SLAMS can be efficient elec- 
tron accelerators dMann & Classed 1 19971) . Their magnetic am- 
plitude reaches several times that of the background field and 
they can propagate wit h a super-Alfvenic speed, because the y 
convect with the ions (iMann et all 11994 iBehlke et alj |2003|) . 



SLAMS have also been observed in simulations (IScholer et al.l 
|2003|) . The acceleration of electrons by sub-structures of quasi- 
parallel shoc ks, which may be S LAMS, has been observed in the 
solar corona dMann et al.lll998l) . 

The absence of self-consistent kinetic models of oblique 
shocks implies that they can curr ently be studied only numeri- 
cally with particle-in-cell (PIC) (Dawson 1983; Eastwood 199 lb 
orwiffi Vlasov simulations ( Arber & Vann 2002; Sircombe e t al] 
l2006l) . The pioneering PIC simulations of plasma slabs, which 
collide with a sp eed of 0.9c and in the pr e sence of an obliqu e 
magnetic field dBessho & Ohsawal 119991: IZindo et al.1 l2005h . 
have evidenced the formation of a shock that accelerated the 
electrons to ultrarelativistic speeds and ampli fied ffie magnetic 
field. A more recent PIC simulation study dDieckmann et al.l 
l2008bh has shown, that the shock formation is triggered by 
an energetic electromagnetic structure (EES). The simulation 
could demonstrate that an approximate equipartition of the ion, 
electron and magnetic energy densities is established. However, 
these simulations could only resolve one spatial dimension due 
to computer constraints, which is not necessarily realistic for 
mildly relativistic collision speeds. 

We perform here a case stud y with initial cond i tions, w hich 
are similar to those employed by iDieckmann et al.l (l2008b ). We 
reduce the collision speed to 0.5c and lower the temperature. 
The plasma cloud representing the subshell is ten times denser 
than the plasma cloud that represents the ambient plasma (inter- 
stellar medium), into which the shell expands. A guiding mag- 
netic field is quasi-parallel to the plasma flow velocity vector 
and it results in an electron gyrofrequency that equals the elec- 
tron plasma frequency of the dense cloud. These bulk plasma pa- 
rameters have be en selected with the intention to enforce a pla- 
nar shock front dHededal & Nishikawal 120051; iBret etaLl 120061; 
IDieckmann et al.l l2008b'). by which we can model the shock in 
one- or two-dimensions in space. The reduced ion-to-electron 
mass ratio we use allows us to model this collision in form of 
a 2(1/2)D particle-in-cell (PIC) simulation, which resolves two 
spatial and three momentum dimensions. 

Our results are summarized as follows: We find higher- 
dimensional structures (density filaments), which initially form 
at the front of the tenuous plasma cloud and expand in time. The 
front of the dense cloud, which turns out to be the most rele- 
vant structure, remains planar and its filamentation is delayed 
but not suppressed by the guiding magnetic field and the high 
temperature. An EES grows ahead of the dense plasma cloud 
before it has become filamentary and the magnetic amplitude of 
the EES reaches a value several times the one of the initial guid- 
ing magnetic field. The EES is pushed by the dense cloud and its 
phase speed in the rest frame of the tenuous cloud is compara- 
ble to the cloud collision speed or twice the Alfven speed in the 
tenuous cloud. The front of the EES expands at an even higher 
speed. Its high speed and amplitude may imply that the EES is 
a SLAMS. The front of the dense cloud and, consequently, the 
EES are slowed down by the electromagnetic wave-particle in- 
teraction. It is thus not possible to define a rest frame moving 
with a constant speed, in which the EES is stationary. This would 
be necessary to measure the frequency of the EES accurately. 
However, the amplitude distribution of the EES suggests that 
its oscillation frequency is below the ion cyclotron frequency. 
Electromagnetic waves are destabilized by the EES ahead of the 
cloud ove rlap layer thro ugh what we think is a four-wave in- 
teraction dGoldsteinlll978l) . The electrons are accelerated in the 
combined wave fields of these waves and in the forming shock to 
highly relativistic speeds. The simulation shows though that the 
strongest electron acceleration occurs at the position, where the 
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B-Field 
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V, -Vb 

Fig. 1. Initial conditions. Two clouds are modelled. Each cloud 
occupies one half space and their mean velocity vectors v/, and 
-V/, point along the x-direction. The background magnetic field 
B() II (1, 0, 0.1). A convection electric field points along ^y. 



shock-reflected ion beam is forming. The plasma collision re- 
sults in a substantial electron acceleration and also in an amplifi- 
cation of the magnetic energy density by one order of magnitude 
within the EES and the forming shock. Both values are probably 
limited by the reduced ion mass of 400 electron masses, which 
we must employ. Radiative processes, which are not resolved 
by the PI C code, will at this stage star t to in fluence the shock 
evolution (ISchlickeiser & Lerchd 20071 l2008l) and we stop the 
simulation. 

This paper is structured as follows. Section 2 discusses the 
particle-in-cell simulation method and the initial conditions. 
Section 3 presents our results, which are discussed in more detail 
in Section 4. 



2. Initial conditions and simulation method 

We model the collision of two plasma clouds. Each cloud con- 
sists of co-moving and equally dense electron and ion species. 
This system is sketched in Fig. [1] The dense species represents 
in this model a fast subshell of the SNR blast shell and the tenu- 
ous species the ambient medium. 

The number densities of the electrons with the charge -e and 
mass Me of the dense cloud 1 and of the tenuous cloud 2 are Ugi 
and ne2 - n^i/lO. The number densities of the ions with the 
charge and mass e and m, of cloud 1 and 2 are «,i = n^x and 
na - ne2- The cloud 1 occupies the half space x < xq and the 
cloud 2 occupies x > xq. Both clouds collide at the initial contact 
boundary, which we set to jcq = 0. All four species are initial- 
ized in their respective rest frame with a relativistic Maxwellian 
distribution that has the temperature T -25 keV. This value ex- 
ceeds by far the temperature of the ambient medium (ISM) into 
which the blast shell expands, which is typically less than leV. 
The upstream temperature may, however, be increased by the in- 
teraction of the ISM with the cosmic ray precursor of the main 
blast shell. The temperature T is also higher than that down- 
stream of the SNR shock, which is usually a few keV. The pur- 
pose of selecting this high temperature of the simulation plasma 
is to improve the numerical efficiency. It maximizes the grid cell 
size, which must be comparable to the plasma Debye length, and 
the time step that is connected to the cell size, while it ensures 
that no highly relativistic particles are present to start with. 



We set Vb - 0.268c and the cloud collision speed Vc - 
2vi/(l +v\lc^) is Vc — cjl. The thermal speed vvi = {ksT Ime)^^'^ 
and Ve2 = Vei of the electrons of the clouds 1 and 2 is v^i 

0. 825 Vb- The jump of the electron's mean velocity at x - xq is 
thus comparable to the thermal speed and the discontinuity in 
the electron phase space distribution is not strong. The thermal 
speeds v,i = vq of both ion species are v/i = {nielmif'^Vfi. The 
guiding magnetic field Bq - o(l, 0, 0.1) with eB^o/'Wc - ^ei, 

where D.e\ - (neie^/OT^eo)'^^ is the electron plasma frequency 
of cloud 1 . The convection electric field with the modulus Ec - 
VbB 7,0 with Bj o = B^fl/lO points along y and changes its sign 
at X = xo. The Alfven speed in cloud 2 is ^ 0.23 in units 
of c and the collision has an initial Alfven Mach number ^ 2. 
The yS = 2jjQ(nei + nj\)kBT jB^ - 0.2, where T is given in units 
of Kelvin, and the plasma is strongly magnetized. However, in 
the box frame of reference, the average kinetic energy density 
of the plasma exceeds the magnetic energy density by the factor 
8.4 and the kinetic energy thus still dominates. 

Typical particle number densities of the ambient inter- 
stellar medium, into which the plasma s ubshell would ex- 
pand, are 0.1 cm~^ < no < 1.0 cm dAcero et al ] I2007t 

lEHison & Vladimirovll2008l) . This ambient medium would cor- 
respond to the cloud 2. These densities are not always known, 
they change as a function of space and are not the same for all 
SNRs. We will thus use normalized units that permit us to scale 
the simulation results to any no- 
Variables in physical units are denoted by the subscript 
p. The time and space are scaled to f = ^e\tp and x - 
/l7'jCp, where As - c/Q^i is the electron skin depth of cloud 

1. The ion skin depth equals 20/lj for the considered ion-to- 
electron mass ratio. The electric and magnetic fields are scaled to 
Ep(AsX) = Q.t,\cmt,e^^ E(x) and Bj,{AsX) = Q.g[mee^^ B(x). The 
charge, the particle number density and the current are normal- 
ized as ppiAsX) = erifipix), np(AsX) = nein(x) and J p(AsX) = 
erieicJix). The normalized Maxwell's equations are 



dB 

dt 



VxE^-— , V-B = 0, 



(1) 



(2) 



VxB^^+J , W-E^p. 
at 

A PIC code solves the Klimontovich-Dupree equations dPupred 
1963) that are derived from the Vlasov equation with the help of 
the method of characteristics. The equations of motion 



^^^qr[E{Xj)+VjXB{Xj)), 



Axj 

PJ = m,.r(vj)vj, — = Vj 



(3) 
(4) 



are solved for an ensemble of computational particles (CPs), 
where the subscript j denotes one CP of the species r. The charge 
q,- and mass of the CP are given in units of e and m^. and 
the momentum is normalized to c. The macroscopic current / 
is obtained from interpolating the current contributions of each 
CP to the grid and from the summation of these interpolating 
currents over all particles. The macroscopic / is used to evolve 
E and B one step in time. The momentum pj of each CP with 
index j is then updated with these new electromagnetic fields, 
which are interpolated to the position jc, of the respec tive CP. 
Interpolation scheme s are discussed b y iDawso i3 (ll983h and the 
one our code uses bv'Eastwoodl (ll991h . 

Our code is relativistic and electromagnetic and we use pe- 
riodic boundary conditions in all directions. The simulation rep- 
resents the x,y plane. The simulation box length is - 5330 
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Fig. 2. Electromagnetic field amplitudes in the simulation at the time T\ . The upper row shows the magnetic field components 
Bx(.x,y) (a), By{x,y) (b) and (c). The lower row displays the Ej^(x,y) (d), Ey{x,y) (e) and E^(x,y) if). 



along X, which is resolved by 1.8 x 10"* grid cells. The length 
Ly - 60 along y is resolved by 200 grid cells. Cloud 1 occu- 
pies the x-interval -Li72 < x < 0, while cloud 2 occupies 
< X < L J2. No new particles are introduced at the bound- 
aries. The clouds convect away from the initial contact boundary 
six — with the speed modulus vi, and the simulation is stopped 
before the rear boundaries of the clouds enter the area of interest. 

Each plasma species is represented by 160 CPs per cell. The 
density difference of both clouds thus implies a lower weight 
for the CPs of cloud 2. We employ a mass ratio mj/mc = 400 
and let the simulation run for the duration T,,,,, - 3500. If no 
strong instabilities develop, the particles could move on average 
a distance vi,T sm — 950. 



3. Simulation results 

In what follows, we analyse the data at two times. The field and 
particle distribution is considered at the early time T\ - 1150. 
We examine, which conditions eventually result in the growth of 
the shock. The shock is forming at the later time r^,,,, and we 
determine the plasma state that is reached by the release of the 
ion kinetic energy. 

3.1. Early simulation time 

Figure |2] displays the fields at the time T\ . The spatial interval 
-300 < X < 300, in which the clouds overlap by vi,T\ ^ 300, 
shows filamentary structures in Bx, B, and . We discuss these 
first. The filamentation instability separates the currents due to 
both interpenetrating clouds and a current system develops in 
the simulation pl ane. It drives th e B , which gives rise to the 

through linear ( Tzoufras et aL,20 06) and nonlinear processes 
jDieckmann et al.l 20091) . The B, is modulated along y with a 
wavelength between 5 As and 20/ls, the latter being equal to the 
ion skin depth. At least the largest structures of B. seem to be 
tied to ion filaments. 

Initially the convection electric field along y ensures that the 
magnetic field is co-moving with the plasma. The gyrocenters of 
the plasma particles do not move relative to Bq and no net cur- 
rent is present. Once the plasma clouds overlap, there will be a 
net particle motion relative to Bo and the plasma particles can 
be deflected into the y, z-plane. The derivative along z vanishes 
in our 2D PIC simulation. The Bx can only be generated through 
dxBy - dyBx - d,E, + Jy, which causes structure formation in 
Bx and By. Figure |2j a) reveals that the filamentation of Bx due 
to such currents is strongest close to the front of cloud 2, that 
is for X < 0. The plots of By and E^ do not yet show a clear 
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Fig. 3. Comparison of the planar electromagnetic fields along x 
for y = 0. (a) compares cB, with cBy in units of cBxS), while (b) 
compares cB^ with \Ey in the same normalization. 



filamentation but one has to keep in mind that the apparent rela- 
tive weakness of the modulations of By compared to those of Bx 
is influenced by the different colour scales. The fluctuations of 
Ex in the interval -300 < x < 300 have no visible counterpart 
in By and B-^. The polarization of the electric field fluctuations 
along Vi suggests a charge density wave, but their correlation 
with the small-scale fluctuations of Bx for -300 < x < -100 
demonstrates that they are not purely electrostatic. 

The dominant fields of the energetic electromagnetic struc- 
ture (EES) are observed in Fig. |2] ahead of the dense cloud 1 
for X > 300. The EES is planar. The field amplitudes along y, z 
are practically constant along y. The By and B^ in Figs. |2lb,c) 
reach peak amplitudes exceeding that of the guiding magnetic 
field in this interval. The Ex and Bx at x ^ 300 show an oblique 
modulation. 

We exploit the planarity and constancy as a function of y 
of the amplitudes of By, B^ and Ey for x > 250, where we find 
the EES. Figure |3] plots these amplitudes as a function of x for 
y - 0. The considered field components have their peak values 
at the front of cloud 1 located at x x 300. The amplitudes de- 
crease exponentially for x < 300 and By and B~ reveal a phase 
shift of 180°. The amplitudes of By, B~ and Ey decrease also 
for increasing x > 300, although more slowly and the fields os- 
cillate. The phases of By and B. are shifted by w 90° in this 
interval and the EES is predominantly circularly polarized. The 
electron gyrofrequency is comparable to the plasma frequency 
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Fig. 4. Ion densities at the time Ti normalized to n,i . The ions of Fig. 5. Electron densities at the time Ti normalized to i . The 



cloud 1 are shown in (a) and those of cloud 2 in (b). 



electrons of cloud 1 are shown in (a) and those of cloud 2 in (b). 



for X < 300 and significantly larger for x > 300. The jump in the 
ratio of the characteristic resonance frequencies at x w 300 will 
alter the wave dispersive properties of the plasma, explaining the 
qualitative change the fields of the EES undergo as we cross this 
boundary. 

The By and of the EES have a left-hand polarization for 
X > 300. Consider a coordinate system defined by x, By and Br. 
The By has a maximum at x « 340, where B- = 0. As we go to 
increasing values of x, the amplitude of By decreases and that of 
B^ becomes negative. The magnetic field vector rotates counter- 
clockwise. A deviation from a circular polarization is provided 
by B. Q + and by + 0, but their amplitudes are relatively 
low compared to B, and B, for 300 < x < 600. 

The B, is almost in phase with Ey. The comparison of B, with 
and the comparison of B- with for x > 300 would show 
that they are shifted by 90°, while By is shifted by 180° rela- 
tive to E,. Figure|2]shows this most clearly at x x 500. The wave 
components = (0, Ey, E,) and = (0, B,,, B,) perpendicu- 
lar to Vh can thus be connected through x -Vj, x Bj^, where 
we assume that Vb II x. The EES and its magnetic field is thus 
convecting with the plasma of cloud 1 . The wavelength of E^ in 
Fig. |2] along x equals that of all field components except B ^ for 
X > 300. Since all components of E are strong, the EES has an 
electric field component aligned with Bq- The EES has the wave- 
length Aw ~ 120/t.s, coinciding with x Inimilnief^^ As- This 
Aw, together with a convection of the EES with the Vb of cloud 
1, points at a connection of the EES to phase space structures in 
the io n distribution, which is typical for SLAMS (iBehlke et alj 
120031) . 

The magnetic oscillations observed in Fig. |2] in the spatial 
interval -300 < x < 300, in which both clouds overlap, must be 
tied to currents. The currents originate from a spatial redistribu- 
tion of the plasma cloud particles and the latter can be examined 
with the help of the density distributions of the four plasma pop- 
ulations. 

Figure |4] displays the densities of the ions of cloud 1 and 
of cloud 2. The ion density of cloud 1 has a peak at x » 300. 
This high-density structure is planar and shows no signs of fila- 
mentation along y. The filamented front of cloud 2 at x -300 
shows no high-density structure. The ion densities of both clouds 
are only strongly modulated well behind the expansion front of 



cloud 1 . Modulations of the density of cloud 1 of up to 40% are 
observed for -200 < x < and the density of cloud 2 is mod- 
ulated even stronger for -300 < x < 0. Weaker modulations of 
both ion clouds extend up to x 200 and their wavelength de- 
creases with increasing x. The filaments of both ion clouds are 
aligned with the x-axis in the interval -250 < x < -100. The 
ion density of cloud 1 is highest at y and lowest at y ^ 10, 
while the density of cloud 2 shows the opposite behaviour. The 
ions of both clouds and their currents are spatially separated, as 
expected. The ion filamentation takes place on the same scale as 
the oscillations of B. in Fig.|2|c) in the same interval, reinforcing 
their connection. The ion density modulation of cloud 2 oblique 
to the X-axis in the interval -300 < x < -250 is not balanced by 
the ions of cloud 1 . 

It is remarkable that the amplitude of B. in the interval 
-300 < X < 0, which displays the strongest ion filamentation, 
is an order of magnitude lower than that of the EES. The latter 
must thus have a different origin. 

Figure |5] shows the density distributions of the electrons of 
cloud 1 and 2. The electrons of cloud 1 accumulate at x 300 
and their density is modulated on scales of an ion skin depth in 
the interval -300 < x < 0. Similar structures are revealed by the 
electron density of cloud 2 in this interval. These patterns follow 
closely those of the ions in Fig.|4] which is at least partially en- 
forced by the need to satisfy the quasi-neutrality of the plasma. 
This is evidenced, for example, by the electrons of cloud 1 for 
-300 < X < -250. Their density shows the same oblique strip- 
ing as the ion density of cloud 2 in Fig.|4jb). The density max- 
ima of the ions of cloud 2 and of the electrons of cloud 1, as well 
as their minima, follow each other closely in this interval. It is 
the electrons, rather than the ions of cloud 1, which compensate 
the net ion charge in this interval. Another consequence of the 
quasi-neutrality is seen at the front of cloud 1, where the elec- 
trons reach a high density. This density peak is well-defined and 
confined to a naiTow x-interval, despite the high thermal electron 
speed, which should rapidly disperse them. The high ion density 
in this interval provides the confinement. In contrast, the cloud 
2 shows no ion density enhancement in Fig. |4tb) at its front at 
X ^ -300. Consequently, the electrons of cloud 2 are spread out 
along X and the fastest electrons have reached an x < -500. 
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Fig. 6. Electron densities at the time Ti close to the front of cloud 
1 and normalized to n^i. The electrons of cloud 1 are shown in 
(a) and those of cloud 2 in (b). 
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ties are multiplied by the factor 4. 



The ions and electrons of cloud 2 in Fig.|3b) and Fig.|5fb) 
show patterns in their density ahead of cloud 1 for x > 300. 
These patterns resemble those in and at x w 300 in Fig.|2] 
The box width along y upstream of cloud 1 is approximately one 
ion skin depth of cloud 2, because rtn/na - 10. The structures 
just ahead of x = 300 perform two oscillations along y. These are 
potentially filaments, but they cannot be caused by the interac- 
tion of cloud 1 and 2 because they are outside the cloud overlap 
layer. The spatial modulation provides a hint. The filament start- 
ing at X ^ 300 and y 30 in Fig. |5lb) changes the value of y, 
at which the density is highest, continuously as a function of x. 
One oscillation takes place on a distance » 120 along x, which 
is the wavelength Aw of the EES. 

The electrons of cloud 2 in Fig. |5lb) close to the front of the 
cloud 1 at X 300 are also modulated on a smaller scale, which 
is akeady evident from Fig.|5] Figure |6] concentrates on this in- 
terval. Filaments can be seen with a thickness ranging from As to 
about 5As close to x 280 and they have vanished already a few 
tens of As behind the front of cloud 1 . Their size and even more 
so their short lifetime due to a rapid thermalization indicates that 
these structures are caused by a filamentation of the electrons. 
Figure HJb) confirms this, because the large ion filaments form 
far behind the front of cloud 1 , implying that the ions have not 
been thermalized by their passage through the interval depicted 
by Fig.|6] 

The densities of the four species in Fig.|2]are integrated over 
y. The plasma density of cloud 1 reaches 1.5 times its initial 
value at its front at x x 300 and it rapidly vanishes for further 
increasing x. No ions of cloud 1 and almost none of its electrons 
have thus propagated upstream of its front. The plasma density 
of cloud 2 is almost constant for x > 300. The electron gyrora- 
dius in is less than As and they should fall behind the ions 
of cloud 2 once they encounter the EES. One may expect that 
the electrons of cloud 2 accumulate at this position, this is how- 
ever not seen in Fig. [T] The absence of electrons of cloud 1 at 
this position and the need to satisfy the plasma quasi-neutrality 
transports the electrons of cloud 2 across B^, which will accel- 
erate these electrons orthogonally to x. 

The ion density at the front of cloud 2 and the electron den- 
sity of cloud 1 at X « -300 are weakly enhanced. These in- 



creased densities are caused by the oblique structures in the in- 
terval -300 <x < -250 in Fig.gtb) and Fig.Ea). The velocity 
vectors of the ions moving in the oblique filaments have been 
deflected away from Vb- This decreases their flow speed along x 
and the ions accumulate. 

The electron densities of both clouds are nonuniform within 
-600 < X < 300, evidencing their interactions through insta- 
bilities. The density dip of the electrons of cloud 2 for 250 < 
X < 300 coincides with the filamentary structures in the same 
interval in Fig.|6] further evidencing that this dip is caused by an 
electron filamentation instability without a significant involve- 
ment of the ions. A notable fraction of the electrons of cloud 2 
have convected to x ^ -600 and some even beyond that. The 
densities of both ion clouds do not yet show a significant modu- 
lation, apart from the accumulation at x w 300. The ions in the 
layer within -300 < x < in Fig.|4]are only redistributed in the 
X, y-plane, but they do not accumulate along x through a shock 
compression. 

The interaction of the particles with the electromagnetic 
fields implies that their dynamics involves at least one spatial 
and three momentum dimensions. We display the phase space 
projections, which we consider to be the most relevant ones. The 
ion phase space density /j(x, p^) can be used to determine if and 
when the shock forms. The ions of cloud 1 and 2 must merge to 
form the downstream region. The electron distribution /.(x, y) 
shows the range of Lorentz factors attained by the electrons. 

Figure [8] displays these distributions. The ions of cloud 1 
convect practically freely. Only the tip of this cloud at x w 300 
is modified. The electron gyrofrequency coce = e\Bo\/mg equals 
Qei, with which we have normalized the time. The ions have gy- 
rated ri/400 w TT or one half of a gyro-orbit in Bq until now. 
The uniform Bq is, however, almost parallel to x and it cannot 
cause by itself the ion rotation at the front of cloud 1 . The mag- 
netic field - (Bj + Bj)^'^ of the EES exceeds B.i_o at the front 
of cloud 1 . The ions will have gyrated by a significant angle in 
this B^, because the amplitude of the EES grows rapidly to a 
large amplitude. The EES also accelerates the electrons of cloud 
2 in the interval 300 < x < 500, which are transported across 
B± by the ions of cloud 2. The ions of both clouds must provide 
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Fig. 8. Pliase space densities at the time Ti . (a) corresponds to 
electrons and (b) to ions. The greyscale denotes the base 10 log- 
arithm of the density. 



through the EES the energy required to accelerate the electrons. 
The modification of the ion phase space distribution can thus be 
attributed to the EES and to the electron acceleration. 

Weak oscillations of the mean momentum {px){x) can be 
seen in the ions of cloud 2 in the interval 300 < x < 600 in 
Fig. [H We can interpret them as follows. The ions of cloud 2 
upstream of cloud 1 move through the EES at the speed ~ -v^ 
in the box frame towards decreasing values of x. A strong cir- 
cularly polarized magnetowave with a wavevector aligned with 
X would force the ions onto a circular path in the Py,Pz plane, 
in which the mean ion beam momenta (pyXx) and {pz)(x) vary 
as a function of x with the wave's periodicity. Since B._o is not 
negligible, this motion must also modulate (p^Xx). Finally we 
see that the ions of the tenuous cloud 2 have been heated up by 
their filamentation at the front of cloud 2 at x w -300. 

3.2. Late simulation time 

Unless stated otherwise, we analyse now the simulation data 
at the time r,,,,, - 3500, when we stop the simulation. If both 
beams had streamed freely, they should have propagated on av- 
erage the distance VhTsm ~ 950. 

Figure |9] displays the densities of both ion clouds. The front 
of the ions of cloud 1 is located at x ^ 750, which is well be- 
hind VbTsim ~ 950. The ions of cloud 1 have thus been slowed 
down and the peak density reaches a value ^ 6 ?A x ^ 750. This 
high-density structure is still planar, but its density is not con- 
stant as a function of y any more. It decreases to 4 at x ^ 750 
and at y Si 40. The front of the ions of the cloud 2 is located 
at X X! -950, which equals the distance expected from the free 
streaming of the ions. Both ion clouds reveal density modula- 
tions on a scale less or comparable to the ion skin depth. The 
largest is found in the ions of cloud 2 at x -300 and the inter- 
val with a reduced density ranges from y to y 30. The ions 
of cloud 2 also show structures in the density distribution for 
X > 750, which is outside the cloud overlap layer The wavevec- 
tor of these oscillations is parallel to x. 

Movie 1 shows the time evolution of the ion densities un- 
til the time Tj,,,,. We first discuss the filamentation in the cloud 
overlap layer The filament formation is seen clearly, when each 
cloud has expanded by the distance 100 or for a time ri/3. The 




(a) X-position: Ions cloud 1 




(b) X-position: Ions cloud 2 

Fig. 9. (Colour online) Ion density distributions at the time T^m 
in units of tin. The ions of cloud 1 are shown in (a) and the ions 
of cloud 2 in (b). 



filaments are asymmetric, a consequence of the unequal cloud 
densities. The ions of cloud 2 show at this time and in the in- 
terval -100 < X < several filaments aligned with the x-axis. 
The filaments have a constant density along y and a width 10. 
These filaments are separated by intervals with a reduced den- 
sity, which have a width of only a few As. The latter have coun- 
terparts in the ions of cloud 1 . The wave vector of these filaments 
is parallel to the y-axis at this time. 

We also observe density fluctuations of the ions of cloud 2 
with a wave vector along x, which develop first in the section 
of the cloud overlap layer with x < 0. This is an instability de- 
veloping as the ions of cloud 2 propa gate through the plasma of 
cloud 1 . A Buneman-type instability teunemanlll958h develops 
between the ions of cloud 2, which can be considered to form an 
unmagnetized beam on electron time scales, and the magnetized 
electrons of cloud 1 . The phase space density oscillations of the 
ions of cloud 2 are periodic with a wave length A % Anvhl^ei 
(not shown), which is typical for the Buneman instability involv- 
ing hot electrons. The phase speed of the unstable wave is close 
to the speed of the ion beam, here composed of the ions of cloud 
2, turning it into a slow oscillation in its rest frame. The ions 
can react to it and they will form phase space holes after a suf- 
ficiently long time. The ions of cloud 1 see the Buneman wave 
as a rapid oscillation to which they cannot react. This explains 
why the ions of cloud 2 in movie 1 at the time T\ and the elec- 
trons in Fig. |5] show the density fluctuations, but not the ions of 
cloud 1 in movie 1. These charge density fluctuations give the 
oscillations of along x in Fig. |2jd). The electrostatic fields 
probably couple to B^, because the Buneman instability is not 
purely electrostatic in the presence of the strong oblique Bq. 

The ion filaments are observed to move in time in the pos- 
itive y-direction, when the clouds have expanded by the dis- 
tance 100 along X. The flow pattern becomes complicated, when 
cloud 1 has expanded to x ^ 400. The filaments in the interval 
< X < 300 still convect to increasing y, while the ion filaments 
of cloud 2 with x < move in the opposite direction. We may 
interpret this observation as follows. Only the ions in the cloud 
overlap layer are not co-moving with Bq, because we cannot de- 
fine through the convection electric field a reference frame, in 
which the ions of both clouds are at rest. The ions in the cloud 
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Fig. 10. (Colour online) Electron density distributions at the time 
Tsim in units of n^i. (a) displays the electrons of cloud 1 and (b) 
the electrons of cloud 2. 



overlap layer rotate with a velocity amplitude set by their mean 
speed relative to the rest frame of the overlap layer, rather than 
by their smaller thermal speed outside it. 

The ions gyrate in the cloud overlap layer orthogonally to 
Bo in a plane that is almost parallel to the py, p, plane. The ions 
of cloud 1 and cloud 2 are counter-propagating and the mag- 
netic field should deflect them into opposite directions. The B, 
separates the filaments of both ion clouds. This implies in a 2D 
geometry that the fila ments are locked as l ong as the magnetic 
trapping force by (^ Davidson et al.|[T972 ') is stronger than the 
Lorentz force due to Bq. A common direction of motion is estab- 
lished. The direction of motion that varies with x and later also 
with y is probably tied to the ion gyro-phase. Movie 1 suggests 
at late times that the filaments of cloud 2 close to its front merge 
and separate again. The B. is here insufficiently strong to sepa- 
rate the ion filaments in the x, y plane. The Bq introduces a 
complicated flow pattern in the 5D phase space spanned by x, y 
and by the three components of p. The filaments do probably not 
merge in this high-dimensional space but they appear to do so in 
the simulation plane. 

Now we turn towards the filaments upstream of cloud 1 . As 
we approach the time Ti, when the clouds have expanded by the 
distance 300 along x in the movie 1, the filaments in the cloud 
overlap region have fully developed and we observe first signs 
of a filamentation of the cloud 2 ahead of the front of cloud 1 . 
These filaments are initially oblique as Fig.|4jb) depicts. When 
cloud 1 has reached x x 600, the upstream region ahead of cloud 
1 shows strong density oscillations of the ions of cloud 2 along 
the X-axis, which are also depicted in Fig.|9|b). They cannot be 
related directly to the EES, because their wavelength along x 
is about Aw/3. The patterns and their time-evolution eventually 
become complicated and presumably also high-dimensional. 

The density distributions of the electrons of both clouds are 
correlated with those of the ions in the overlap region -950 < 
X < 750, as the Fig. [10] is evidencing. An example is here the 
density pattern, which crosses y = at |x| < 100 and is visible 
in all four species. The filaments in the interval -500 < x < 500 
are, however, more diffuse than those of the ions in Fig.|9l The 
density peak of the electrons of cloud 1 at x » 750 is comparable 
to that of the ions. Again, the electrons are confined to maintain 
the quasi-neutrality of the plasma. The front end of the electrons 



of cloud 2 has been replaced by a gradual decrease of the elec- 
tron density between -1500 < x < -800. The electrons with 
X < -950 are the fastest ones of cloud 2, which have outrun the 
bulk of the electrons moving with -Vh- The ripples in the elec- 
tron density of cloud 2 at .jc > 750 are practically identical to 
those of the ions in Fig.|9fb). 

Movie 2 reveals the time-evolution of the electron densities. 
Density structures form rapidly close to the initial cloud contact 
boundary. These structures have no counterpart in the ion distri- 
butions in movie 1, until the cloud overlap layer covers |xi $ 50. 
This initial plasma evolution is thus determined by the electron 
filamentation. The density of the electrons of cloud 1 is enhanced 
in the cloud overlap layer, while the electron density of cloud 2 
is decreased (See also Fig.|2ta)). Already at this early time, the 
fastest electrons of cloud 2 outrun the bulk of the electrons to 
form a low-density tail along x and one consequence is the re- 
duction of the electron density of cloud 2 in the cloud overlap 
layer. A strong correlation between the ions of both clouds and 
the electrons of cloud 1 in the overlap layer is observed, when 
the latter covers \x\ ^ 100, demonstrating that these three species 
are involved in the formation of the large scale pattern. The elec- 
trons of cloud 2 show a distinct behaviour This is also observed 
at a later time, when the overlap layer covers |x| ^ 200. The elec- 
trons of cloud 2 are redistributed at this time immediately after 
they have crossed the front of cloud I at x ^ 200, which is also 
shown by Fig.|6l and their density is not visibly correlated with 
those of the other species. 

This behaviour can be understood qualitatively. Consider a 
dense electron beam and a tenuous electron beam that undergo 
the filamentation instability, which separates their currents. A 
highly nonlinear filamentary electron distribution of the tenuous 
beam can be compensated by a weak perturbation of the dense 
beam. We thus expect that the spatial density modulations of the 
tenuous beam become nonlinear first. We also expect, that the 
high ion inertia results in a delay of their filamentation compared 
to that of the electrons and both species of cloud 2 will interact 
differently with the plasma of cloud 1 . 

The ions of cloud 2 in Movie 1 show structures similar to 
those in cloud 1, which have initially no counterpart in the elec- 
trons of cloud 2 (Movie 2). The ions of cloud 2 reveal a higher 
degree of non-uniformity than the electrons of cloud 1, when 
the overlap layer spans the interval \x\ ^ 200. We can attribute 
this to their lower mean density. The increasingly pronounced 
large-scale modulation of the plasma density of cloud 1 at the 
time Ti imprints itself onto the electrons of cloud 2. The latter 
show a modulation on two scales in Fig. |5] The density evolu- 
tion of the electrons and of the ions eventually becomes similar, 
except close to the front of cloud 1 . Far behind this front, the 
electrons have thermalized. Since they have no free energy left 
which could lead to independent behaviour their density distri- 
bution follows that of the ions. 

We explore now the density evolution of the electrons and 
ions of cloud 2 ahead of the front of cloud 1 in the movies 1 
and 2. We determine now, if we find multiple particle beams 
upstream of cloud 1 that could give rise to beam-driven instabil- 
ities. The electrons of cloud 2 thermalize after they have crossed 
the front of cloud 1 in Fig.|6]and we can no longer distinguish be- 
tween the electrons of either cloud. Movie 2 shows no electrons 
of cloud 1 leaking upstream and the same should hold for the 
electrons of cloud 2, which have been mixed with those of cloud 
1. This is supported by Fig.|2l which demonstrates firstly that no 
significant numbers of electrons from cloud 1 diffuse upstream 
and, secondly, that the electrons of cloud 2 ahead of cloud 1 
have a uniform density and convect with the ions to preserve the 
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Fig. 11. Electromagnetic field amplitudes at the time Tsim- The upper row shows the magnetic field components (a). By (b) and 
5^ (c). The lower row displays the (d), Ey (e) and E^ (f). 



quasi-neutrality of the upstream plasma. Figure |8] demonstrates 
that no ions of cloud 1 leak upstream of its front and that no ions 
of cloud 2 are reflected at x = 300 to form a beam upstream of 
cloud 1. Movies 1 and 2 together with Fig.|2]evidence, however, 
the growth of the EES and the onset of the filamentation up- 
stream of cloud 1 prior to the time T[ . These two processes thus 
grow in an interval, in which we only find the plasma of cloud 
2 and no further particle beams. We may thus exclude for our 
case study and until the time Ti that particle beams (which can 
be formed by leaking downstream plasma or reflected fast ion 
beams in a quasi- parallel magnetic field geometry or in unmag- 
netized plasmas ( Malkovlll998t iMartins et aLll2009l) ). result in 
the growth of the EES and in the observed upstream instability. 

Figure [TT]displays the fields at the time Tsim- The nonplanar 
B^ reaches amplitudes comparable to B.v.o/2. Both components 
of evidence planar oscillations with an amplitude exceeding 
Z? t and with a wavelength » 120. The wavelength is unchanged 
and the EES is still circularly polarized. The amplitudes of By 
and B~ peak in the interval 600 < x < 1000 close to the front of 
cloud 1. Evanescent waves reach far upstream up to x ^ 1700. 
The Ej_ is more diffuse and Ey shows structuring along y. Quasi- 
planar structures close to and ahead of the front of cloud 1 are 
revealed by E^ for 700 < x < 1000. The wave vectors of these 
oscillations point predominantly along x and their characteristic 
wavelength is short compared to that of the EES. They resemble 
the structures of cloud 2 upstream of cloud 1 at late times in the 
movies 1 and 2. 

The evolution of the fields is presented in the movie 3 for 
the three components of B and in the movie 4 for the three com- 
ponents of E. Until the time Ti, when the overlap layer covers 
|x| < 300, the fields can be subdivided into those in this layer 
and those of the EES. 

We discuss first the electromagnetic fields in the cloud over- 
lap layer. The filamentation is observed best for B^ and Ey. 
The B^ is the dominant component of the filamentation insta- 
bility in our 2D geometry and for a flow veloc ity vector along x 
dBret et al.ll2008t ISchhckeiser & Shuk la 2003). This component 
indeed shows a filamentation, but the adaptation of the colour 
scale to the EES reduces the contrast and it is not clearly visi- 
ble. The B V, Ey and B- show filamentary structures with a thick- 
ness below or equal to the ion skin depth. The filamentation of 
these field components shows the same flow patterns as the ion 
distribution in movie 1 . The ion beam filamentation is thus the 
primary driver of the fields. Oscillations of E^ with a short wave- 
length along X are observed, which we have related previously 
to a Buneman-type instability. 

The EES starts to grow immediately at the cloud collision 
boundary and it is thus a consequence of our initial conditions. 



The initial jump of the convection electric field £, across x = 0, 
its constant amplitude along y and the vanishing derivative along 
z imply that dE^/dx - -diB.. A seed field for magnetowaves de- 
velops at the initial contact boundary. We may interpret its con- 
sequence as follows. The magnetic field is approximately sta- 
tionary in the dense plasma. This can be seen from the Ey in the 
cloud overlap layer in Movie 4, which is practically identical to 
that in cloud 1 outside the layer. The deflection of the electrons 
of cloud 2 by this B., which is localized in x, will be stronger 
than the deflection of the ions of cloud 2 and 7,. 9^ 0. The EES is 
unstable, which implies that the current amplifies the magnetic 
perturbation. The EES reaches a B± ^ B^o akeady when the 
cloud overlap layer has expanded to cover |x| ^ 100. The EES is 
thus strong enough to cause the deflection of the ions at the tip 
of cloud 1 at the time Ti in Fig. [8] The amplitudes of By and 
grow and the EES expands steadily until secondary instabilities 
occur. 

When the cloud overlap layer covers |x| ^ 500, the struc- 
tures in B^ upstream of cloud 1 in movie 3 start to accelerate 
along y. This acceleration is most easily visible close to the po- 
sition, where the amplitude of the EES peaks. The time, when 
the structures of B^ start to accelerate along y, coincides with 
a similar flow of the structures of Ey and E, in movie 4. Prior 
to their acceleration, the structures in B ^ resemble the oblique 
structures in the plasma density of cloud 2 upstream of cloud 1 
found at the time Ti in the Figs. |4] and |5] The density structures 
and those of By are modulated along x on scales comparable to 
the wavelength of the EES. 

We may interpret the formation and the acceleration of the 
field structures in terms of an acceleration of the incoming par- 
ticles of cloud 2 in the wave field of the EES. If its wave field is 
not perfectly uniform along y, then the particle deflection is not 
equally efficient for all y, which will cause density modulations. 
The particles of cloud 2 with their spatially nonuniform den- 
sity are then accelerated in the y,z plane by the EES, resulting 
in nonuniform currents. A spatially nonuniform current could 
drive the nonplanar field structures, which we observe in movies 
3 and 4. A nonuniform J, would, for example, yield structures 
in By, By and E^ by dxBy - dyBy - d,E, + J,. The motion of B^ 
and also of E^ and By would then be induced by that of the den- 
sity structures. The growth and the accelerating motion along y 
of the structures are presumably tied to the increasing ampUtude 
of the EES. 

Movie 4 shows that the E^ oscillates outside the cloud over- 
lap layer. The amplitude of E^ is constant within each cloud on 
scales much larger that and the oscillations start instantly. 
Thus, they cannot be connected to the processes at x a; 0, be- 
cause of the finite light speed. These oscillations are periodic in 
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Fig. 12. The Bx.{x,y - 0, f) in the reference frame of the simula- 
tion box. The By is shown in (a). The two white lines correspond 
to X = Vi,f and x - 0.8cf. The B, is shown in (b). 
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the electron densities and (b) the ion densities. The solid curve 
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time and their frequency is the electron plasma frequency of the 
respective cloud. These oscillations presumably originate from 
a slight initial mismatch of the currents of electrons and ions 
along X in each cloud, which results through d,Ey: - -J^ in elec- 
trostatic oscillations at the local plasma frequency. 

Structures in the plasma of cloud 2 with a wave vector that 
is aligned with x were observed in the Figs. l9l and [TOl for 750 < 
X < 1200. Movies 3 and 4 provide additional information. These 
short-scale structures do not show up in B^ and hence they do 
not have a strong electromagnetic component. The modulations 
are most clearly visible from for example in Fig. [TTl d). A 
polarization along Vh indicates that these waves are electrostatic. 
The density structures in cloud 2 should thus be charge density 
waves, similar to the Buneman waves we have observed in Fig. 
|2fd). This charge density wave is not the only one upstream of 
cloud 1 . Movie 3 shows at late times and for large x > 0, that 
magnetowaves with a longer wavelength than that of the EES are 
generated that move faster than the EES. 

A wave conversion in form of a four-wave coupling of a 
cir cularly polarized, large amplitude wave has been proposed 
by iGoldsteinI ( Il978l) . The four interacting waves all have par- 
allel wave vectors and they will thus all result in plasma and 
field modulations along x in our geometry. The four-wave in- 
teraction results in a forward-propagating magnetowave with 
the same polarization as the driving wave, which might be the 
long- wavelength magnetowave observed in movie 3. It is faster 
than the super-Alfvenic EES, implying that it should propagate 
in the high-frequency branch of the left-hand polarized mode 
dTreumann & Baumiohannlll997l) . This parametric wave inter- 
action also pumps a charge density wave, which could result in 
the density modulation observed in the movies 1 and 2 ahead of 
cloud 1 . Finally, the interaction should produce another magne- 
towave that moves in the opposite direction than the EES. 

We could not identify this fourth wave in E(x, y, to) and 
B(x,y, to) at any given time to- This either means it is not there 
or that its amplitude is close to noise levels. Its wave vector is 
aligne d with that of the EES and it is a magnetowave (Goldstein 
119781) . It should thus be detectable in slices of Bj_(x) for y - 0, 
which we can plot as a function of time. The space-time corre- 



lation allows us to identify structures, which are not easily seen 
from individual time slices. Figure [12] shows them as a function 
of time. Both components of B^ show at x x 800 and t x 2700 
waves propagating antiparallel to the EES. These fast waves 
originate from processes taking place at the rear end of cloud 
2, which are fuelled by the different expansion speeds of elec- 
trons and ions into a vacuum. These waves can be seen already 
before they start to interact with the EES. The By amplitude is 
modulated at x « 1700 and t x 1 100. These waves are amplified 
as they pass through the EES. We have to emphasize though, 
that the backward-propagating magnetowave is only barely visi- 
ble in this plot and that we have not shown that its amplification 
by the EES results from a coherent 4-wave interaction. The lat- 
ter would require more detailed correlation studies of the wave's 
phases. Our current simulation data is too noisy and we leave this 
to dedicated ID PIC studies that can employ a higher number of 
particles per cell to reduce noise levels. 

Two lines are fitted in Fig. [T2l a). The x - vi,t is approxi- 
mately twice the Alfven speed Va of the cloud 2. Initially, the 
rear end of the EES moves with vi, and it is thus connected to the 
front of cloud 1. The speed of the EES decreases steadily and 
its rear end at the time r„„, coincides with the front of cloud 1 
at X ^ 750 rather than vi,Tsim ~ 950. This deceleration makes 
it difficult to determine the frequency of the EES, which should 
be measured in its rest frame. The first maximum of By moves 
with the front of cloud 1 and shows no change in time. We may 
conclude that its oscillation frequency is low or zero in this de- 
celerating reference frame. The B, shows a change at the front of 
cloud 1, which may evidence a frequency ojees of the EES. 
It may, however also be related to the formation of the shock. We 
can thus only conclude here that loees ~ ^tt/Tsu,, ~ (me/m,)^)^. 
The second fitted line in Fig. [T2? a) corresponds to x = 0.8cf 
and it is an approximate boundary for the front of the EES. This 
expansion speed is 0.87c in the reference of cloud 2. 

Figure [T3] shows the densities of both clouds, averaged over 
y, at the time r,,-,,,. We consider only x > 0. The densities of all 
four species are practically constant for x < 600 and the density 
of each species is comparable to its initial value. The processes 
in the cloud overlap layer have reached a state that does not in- 
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X-position: Ions 

Fig. 14. Ion phase space density at the time Tsim- The grayscale 
shows the 10-logarithmic density normaUzed to the density peak. 



volve a density modulation. No shock is present in this interval, 
because that would result in a plasma compression. The plasma 
density of cloud 1 rises in the interval 600 < x < 750 and de- 
creases to zero for x > 800. The front of cloud 1 is thus still 
well-defined. The density of cloud 1 increases to about 3-4 times 
its initial value in the interval 650 < x < 750, reaching a peak 
density of 6. The plasma density of cloud 2 reveals substantial 
oscillations for x > 800 upstream of cloud 1 . The electron and 
ion density oscillations, which could be observed in the movies 
1 and 2, and the y-averaged densities do not match perfectly. The 
resulting net charge modulation results in charge density waves 
and the £, in Fig.[nid). 

The ion phase space density distribution in the x, plane 
and at the time Tsim, which has been integrated over all other 
position and momentum components, is shown in Fig. [14] The 
ions of both clouds have started to mix in the x, p^-space, which 
is the condition for a shock formation. The ion clouds merge in 
the interval 600 < x < 800, which coincides with the interval 
that has an increased plasma density in the Fig. [T3] The den- 
sity increase is thus due to the shock compression and it will 
result in a downstream region that expands along x. This future 
downstream plasma has not yet thermalized, explaining why the 
density in Fig.[T3]is not constant within 600 < x < 800. 

The forward shock, which will be moving to higher x, will 
be somewhat slower than the collision speed of c/2. It will move 
faster than in the reference frame of cloud 2. This is because 
the plasma has a net momentum in the simulation frame due to 
the higher density of cloud 1 . The shock moving in the direction 
of decreasing x will move slower than vt, in the reference frame 
of cloud 1 . We can already see the slowdown of the front of cloud 
1, which coincides with the amplitude maximum of the EES in 
Fig- El which will form the forward moving shock. 

Figure [14] also reveals oscillations of the mean momentum 
{px){x) of cloud 2 in the interval 800 < x < 1000. The waves that 
can cause such strong oscillations must yield clearly detectable 
field oscillations in the same interval. They must be caused either 
by the charge density waves and the associated E ^ or by the EES 
ahead of cloud 1 , or by the combination of both. The wavelength 
of the oscillations is well below the wavelength » 120 of the 
EES and comparable to that of the oscillations of in Fig.[TTJd) 
with a wavelength A^l'i- The electrostatic potential of the charge 




X-position: Electrons 

Fig. 15. Electron phase space density at the time Tsim- The 
grayscale is the 10-logarithmic density normalized to its peak. 



density waves will indeed cause particle oscillations in the x, p ^ 
plane. The main effect of the EES, namely the gyro-bunching, is 
not necessarily visible from the phase space density distribution 
in the x, plane. The of the EES forces the ions on a gyro- 
orbit orthogonal to x. The electrostatic component of the EES 
provided by E^ will, however, modulate the ion distribution in 
the X, p j-plane. No clear modulation of the ion paths on scales 
~ Aw is visible in Fig. [T4l though. 

The electrons are accelerated to highly relativistic speeds in 
the interval in which we find the charge density waves and the 
forming shock. The peak Lorentz factor y is 120, according to 
the Fig. [15] This peak Lorentz factor is remarkable in that it 
exceeds even the (m,7mf)(y[vc] - 1) » 60, which we would 
obtain if the energy of an electron would equal that of an ion 
moving with the cloud collision speed Vc - cjl. This energy 
is, however, reached only by a tenuous subpopulation of the 
electrons in a limited spatial interval. This strong electron ac- 
celeration is presumably achieved by those fields, that give rise 
to what will become the shock-reflected ion beam. The ions of 
cloud 2 at X «i 750 in Fig. [14] are accelerated from the speed 
-Vh to the speed 0.5c > Vh at this position and their velocity 
thus changes by more than vv. We also observe that the elec- 
trons are not only accelerated in the forming downstream region 
600 < X < 800, but they appear to be pre-accelerated in the in- 
terval 800 < X < 1100. The electron distribution shows narrow 
spikes in this region. The electron phase space density is smooth 
in the interval 500 < x < 750 and the energies are here compa- 
rable to the energy of an ion moving with vv. The peak energy 
of the electrons decreases steadily, as we go from x » 500 to 
X ^ -1500. The characteristic distances crossed by electrons 
with y > 5 will not differ much during the short time Tsim- This 
distance is thus proportional to the time, which has passed since 
they were accelerated. The energy to which electrons have been 
accelerated must therefore have increased with time. 

Figure [16] displays the electron phase space density from 
Fig. [15] which has been integrated over all x to give Niy). It also 
shows the probability N{y) — N{y*)dy* of finding electrons 
with a Lorentz factor y* > y. Most electrons have a y < 5. The 
density of this cool population is exaggerated, because we inte- 
grate over the full simulation box and, thus, over regions which 
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(a) Lorentz factor y 
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Fig. 16. Electron energy distributions, sampled at the time r„„, 
and normalized to their respective peak value, as a function of 
y. The number density N{j) is displayed in (a), while (b) shows 
N{y) = Nif > r) = rN(r)dy\ 



have not yet been affected by the cloud collision. The high en- 
ergy tail in the interval 5 < y < 100 in both distributions can 
be approximated by an exponential and the densities decrease 
rapidly to zero for y > 120. 

Movie 5 provides further insight into the electron accelera- 
tion mechanisms. It is a time-animation of the phase space den- 
sity distributions /i (x, p^) of the ions and fe(x, y) of the electrons 
until the time Tsm- The colour scale denotes the 10-logarithmic 
number of computational particles in units of those of cloud 2. 
One particle of cloud 1 corresponds to 10 of cloud 2. 

During the initial phase of the cloud collision, up to the time 
when the ions of both clouds overlap in the interval S 50 
in movie 5, the ions propagate practically freely. The ion phase 
space distribution at the front of each cloud is increasingly tilted, 
because ions with a higher speed propagate farther in a given 
time interval. The electrons are moderately heated to y » 2.5 at 
the front of cloud 2. The heating mechanisms are the filamen- 
tation instability and the Buneman-type instability, which also 
increase the thermal spread of the ions of cloud 2 at its tip in 
Fig.ISb). 

When the cloud overlap layer covers |x| $ 300, the electron 
heating at the front of cloud 2 has ceased and the electron accel- 
eration sets in at the front of cloud 1 . The slowdown of the ions 
of cloud 1 is visible and it increases steadily, which is consistent 
with a magnetic deflection by the growing EES. When the ions 
of cloud 1 have reached x w 400, the electrons are accelerated to 
y « 6. The movie 5 shows that this electron acceleration is lim- 
ited to the interval covered by the EES. The electrons of cloud 2 
are accelerated as they approach the front of cloud 1 . This accel- 
eration results in a beam in /e(x, y) prior to their encounter with 
the front of cloud 1 and a depletion of nonrelativistic electrons. 
The movie 5 shows this as a detaching of the electron distribu- 
tion from the coordinate axis y = 1 . 

We may interpret this as follows. The electrons of cloud 
2 are dragged with the ions across the Bj_ of the EES and 
all electrons experience the same accelerating force during the 
same time interval, resulting in a similar velocity increment for 
all electrons. The initial electron beam stays compact and the 
electrons move on a corkscrew orbit in the circularly polarized 
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Fig. 17. Phase space distributions at the time O.STj/m. (a) displays 
the electrons and (b) the ions. 



EES jPieckmann et al.ll2008bl) . After this acceleration, the elec- 
trons of cloud 2 cross the front of cloud 1 and they are mixed 
in phase space to form a smooth density distribution. Movie 5 
demonstrates that the electron acceleration by the EES increases 
steadily in time, presumably because their cross-field transport 
across the EES with its increasing amplitude and extent along x 
provides a stronger acceleration for a longer time. 

When the ions of cloud 1 have propagated to x w 500 in 
movie 5, the ions of both clouds start to merge. We show the 
phase space densities fi{x,px) and fe{x,y) of both species di- 
rectly after this merger at the time O.Sr,,-,,, » 2800 in Fig. [17] We 
notice at this time localized and strong oscillations of {px){x) at 
800 < X < 1000 and weaker ones with a shorter wavelength 
within the interval 400 < x < 800. The electrons are accelerated 
to highly relativistic energies by these structures, as expected. 
Consider a force that causes velocity oscillations of the ions of 
Ay =s 0.1c in Fig. [TTl b). The same force will clearly result in 
relativistic electron velocity oscillations. The oscillations of the 
ion mean momentum move in the negative x-direction in movie 
5 and they are amplified as they pass through the EES. 

We may associate these wave structures with the backward 
propagating charge density waves and magnetowaves expected 
from the four-wave interaction, although we have to point out 
that the strong plasma modulations imply a nonlinear stage of 
the plasma, which is not incorporated into the equations describ- 
ing the four-wave interaction. A consequence can be that e.g. the 
strong and compact magnetowave modulates the plasma density 
by the pondero motive force of its fields. The fields resulting in 
the large modulation must belong to a well-defined wave packet. 
Figure [12] evidences at the time 0.8rj„„ the presence of a mag- 
netowave in the interval 800 < x < 1000, in which we find the 
strong momentum oscillation of cloud 2 in Fig.[T7lb). The oscil- 
lations of {px){x) in the interval 600 < x < 700 may correspond 
to the charge density wave, which we observed at late times in 
the movies 1 and 2. Both, the charge density waves and the mag- 
netowave accelerate electrons to highly relativistic speeds. 

We finally assess the energy contained by the fields at the 
time Tsim and the value of plasma collisions of the type we have 
considered here for the magnetic field amplification upstream 
of the main blast shell of a SNR. The magnetic energy den- 
sity in physical units B^/2yUo is normalized to B^/ljUo and the 
electric energy density in physical units eo£'~/2 is normalized to 
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Fig. 18. Field energy densities at the time Tsm- (a) considers the 
magnetic energy density and (b) that of the electric field. 

c^B^/lfiQ. Both are plotted in Fig. [18] Significant field energy 
densities are found only close to the front of cloud I at x ^ 750. 
The magnetic field dominates and its peak energy density ex- 
ceeds that of the electric field by the factor x 25. The slowdown 
of the front of cloud 1 implies that the = "^f^ ^j. decreases, 
where V/ is the speed of the front of cloud 1 that was initially vi,. 
The energy density of the magnetic field is amplified by over an 
order of magnitude by the plasma collision and by the shock for- 
mation. The strongest magnetic fields corresponding to the cir- 
cularly polarized EES coincide with the fastest electrons in the 
simulation. This interval should emit significant electromagnetic 
radiation. 

4. Discussion 

We have described in this paper the collision of two plasma 
clouds at the speed c/2. The ion to electron mass ratio of 400 
has allowed us to model the collision in two spatial dimensions 
until the shock forms. Then we had to stop the simulation. The 
acceleration of electrons to speeds ~ c and the rapid expansion 
speed of the energetic electromagnetic structure (EES) imply, 
that both will quickly reach the boundaries; hence our periodic 
boundary conditions become invalid. Open boundaries would al- 
low the electrons and the wave energy to flow out of the system. 
Howe ver, the instabilities driven by these beams (Martins et al. 
|2009|) and by the EES are important for the upstream dynamics 
and the latter will be adversely affected by open boundaries. 

The filamentary structures can form and merge in the 2D 
geometry we consider here up to the instant when the mag- 
netic repulsion of two filaments with opposit ely directed cur- 
rent enforces their spatial separation (.Davidson et al.l 119721: 
iLee & Lampdl 19731) . at least in the absence of an oblique mag- 
netic field. Then no further mergers occur since only one di- 
mension is available orthogonal to the flow velocity vector. A 
realistic 3D PIC simulation would allow the filaments to move 
aroun d each other and mer ge with other filaments of equal po- 
larity dLee & Lampel[T973l) . A 3D PIC simulation is, however, 
currently impossible for our case study involving ions, because 
of the computational cost involved in resolving ion and electron 
scales. Simulations in three spat ial dimensions are no w feasible 
for the case of leptonic shocks dNishikawa et al.ll2009l) . 
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The densities of the clouds differ by a factor of 10 
and the collision is asymmetr ic, as i n the s imulation by 
Bessho & Ohsawa (1999); Zind oet al.1 ( l2005h ; ISorasio et all 
(120061) . Initially the magnetic field is uniform and quasi-parallel 
to the flow velocity vector Its significant strength together with 
the high plasma temperature of 25 keV and the unequal cloud 
densities reduce the growth rate of the filamentary instabiUties 
(Bret etal.2006,200 7J. 

iDieckmann et al.l (l2008bl) have previously probed the higher 
■y regime appropriate for the internal shocks of gamma-ray 
bursts. Here we consider the mildly relativistic regime, with a 
collision speed 0.5c between both clouds. Such a speed might be 
realistic for a plasma subshell outrunning the main SNR shock. 
Such subshells can reac h relativistic flow spe eds for particularly 
violent SNR explosions dKulkarni et al.l 19981) . In our initial con- 
ditions, the magnetic energy surpasses the thermal energy of the 
dense plasma slab by a factor of 5. The plasma flow implies, 
however, that the box-averaged plasma kinetic energy density 
exceeds the magnetic energy density by an order of magnitude. 
We summarize several aspects of our results. 

4.1. Effects due to initial conditions 

Our initial conditions have resulted in the formation of planar 
wave and plasma structures, the most important one being the 
EES. We think that the EES grew out of a localized seed mag- 
netic field pulse driven by the spatial gradient of the convection 
electric field. The plasma upstream of the dense cloud is desta- 
bilized by this electromagnetic structure and the EES expands at 
the speed 0.87c in the reference frame of the tenuous cloud. The 
energy for its growth and expansion is provided by the kinetic 
energy of the upstream medium, which moves with respect to 
the EES. The shock speed ^ then implies that we have a co- 
herent magnetic layer that expands its width at a speed of at least 
(0.87c - Vc)/(1 - 0.87vc/c) ~ 0.65c, measured in the reference 
frame of the tenuous cloud. It covered about 80 ion skin depths 
at the end of the simulation, showing no signs of a slowdown. 

The EES is a consequence of our initial conditions and 
the growth of the seed magnetic field amplitude could prob- 
ably be reduced but not suppressed by a smoother change of 
the convection electric field, which can be a chieved by a grad- 
ual change of the plasma convection speed (iBessho & Ohsawal 
119991 ; IZindo et al.l |2005|) . However, the seed magnetic field 
could be provided also by waves with a short wavelength, e.g. 
whistlers, and it is thus not unphysical. 

Structures with strong magnetic fields, similar to the EES 
and known as SLAMS, are frequently observed close to quasi- 
parallel shocks in the solar system plasma and they can acceler- 
ate electrons to high energies. They are thus potentially impor- 
tant also for SNR shock physics. Our initial conditions provide 
a possibility to let nonlinear MHD waves grow out of a simple 
simulation setup for a further study. The EES is moving with the 
ions of the dense cloud and it modulates the ions and electrons of 
the tenuous cloud, thereby gaining energy. Its growth probably 
requires an asymmetric plasma collision. 

4.2. Sliock formation 

In this paper we modelled the formation of the shock from the 
initial collision of two plasma clouds. The signatures of the 
shock are evident, including visible thermal broadening behind 
the shock and a dense shock ramp. While filamentation struc- 
tures form ahead of and behind the shock, we note that the 
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Structure is basically planar in the critical foreshock area, where 
electron acceleration is expected to occur This means that one- 
dimensional simulations will be relevant in this region, allowing 
much higher resolution, increased particle number (resulting in 
lower particle noise and a better phase space resolution) and a 
higher ion-electron mass ratio, than is currently found in two 
and three-dimensional simulations. 

It is evident from the simulation that the filamentation is 
not fully suppressed by the guiding magnetic field and by the 
high plasma temperature. Its amplitude has been set such that it 
should suppress the el ectron filament ation if the plasma would 
be spatially uniform ( Bret et al] |2006V This amplitude is appar- 
ently insufficient to suppress the slower filamentation of the ion 
beams and we could even see evidence for an electron beam fil- 
amentation just behind the front. The front of the dense plasma 
cloud maintains its planarity throughout the simulation, but even 
here the density was non-uniform along the boundary. The onset 
of the filamentation was, however, delayed. The likely cause is 
the high density gradient across the front, which alters the elec- 
tron and ion skin depths and thus the characteristic scale of the 
filaments. The gradient is caused by the slowdown of the ions 
by the magnetic field of the EES and by the electron acceler- 
ation. The electrons are confined at the front in the direction of 
the shock normal so that they preserve the quasi-neutrality of the 
plasma, but they can move orthogonally to it. The latter results 
in a drift current. 

4.3. Field amplification 

IVolk et al.1 (12005) and lElUson & Vladimirovl (2008) have de- 
scribed observations of magnetic field amplification above the 
value expected from shock compression in SNRs. At the final 
simulation time the magnetic field energy density is increased in 
strength by over one order of magnitude, exceeding by far that 
expected from the magnetic field compression by the shock. A 
shock compresses only the magnetic field component perpen- 
dicular to the shock normal, which is weak in our case, and the 
amplification of its energy density can reach a factor of 4-7. The 
magnetic energy density at the simulation's end has been com- 
parable to the box-averaged kinetic energy density in an interval 
spanning about 10-20 ion skin depths. The magnetic energy den- 
sity due to the EES was at least twice as high as that of the back- 
ground field in an interval covering 50 ion skin depths. Even if 
we take into account that the kinetic energy density close to the 
shock is increased by the accumulation of plasma, the magnetic 
energy density still constitutes a sizeable fraction of the local 
total plasma kinetic energy density. The EES has provided the 
main contribution to the magnetic energy density and exceeded 
that due to the filaments downstream by two orders of magni- 
tude. 

Throughout this paper, we used normalized quantities in our 
simulation and we can scale the magnetic field amplitude to the 
relevant plasma conditions. If we set the electron density of the 
dense cloud to 1 cm"^, we would obtain a peak magnetic field 
with a strength of 10 mG. However, we have to point out that 
our initial magnetic field amplitude has been higher than that ex- 
pected for the ambient plasma, even if we take into account its 
amplification by cosmic ray-driven instabilities, and our simula- 
tion results may not be directly applicable. 

Where does the extra field come from? Amplification of the 
magnetic field can occur from the electron drift current arising 
from the E x B drift motion in a layer close to the shock that is 
narrower than the ion gyroradius but wider than the electron gy- 
roradius, see, e.g. lBaumjohann & Treumannl(ll996h . The current 



adds to the shock current and increases the jump in the perpen- 
dicular magnetic field. This can only occur when the ion and 
electron gyroradii differ, i.e. not in a pair plasma. The EES has 
a significant /i v-component and \Bj_\ ^ \Bo\. We thus obtain a 
E X B drift orthogonal to the flow velocity vector 

We have also found that the requirement to maintain quasi- 
neutrality of the plasma implies that the upstream electrons are 
dragged with the upstream ions across the EES, which moves 
with the shock. The resulting v x B drift accelerates the elec- 
trons orthogonally to the shock propagation direction, further 
enhancing the net current and the magnetic field. Finally, mag- 
netic fields of SLAMS are provided by the current due to the 
gyro-bunched ions, which rotate in the plane orthogonal to the 
wavevector 

These mechanisms increase the mean magnetic field, and are 
different from the instability described by Bell (2004) who has 
described a cosmic ray streaming instability which can amplify 
turbulent magnetic fields ahead of the shock. We can exclude 
Bell's instability here since we have not found energetic parti- 
cles with a significant density moving upstream, which would 
provide the net current that drives this instability. 

4.4. Electron acceleration and upstream wave spectrum 

The shock retains its planar structure after it forms. A circularly 
polarised large-scale precursor wave, the EES, expands into the 
foreshock. Its wavelength is several times the ion skin depth. It 
gradually rotates the quasi-parallel magnetic field into a quasi- 
perpendicular one at the current layer and it forces the incom- 
ing ions and electrons to interact with it nonlinearly. The ions 
are gyro-bunched and some of the incoming ions of the tenuous 
cloud are reflected by the formin g shock. We hav e found evi- 
dence of a parametric instability (lGoldsteinlll9 78') of the EES 
ahead of the foreshock and we could find at least two waves that 
may be the result of this parametric decay. These waves appear 
at late times, when the EES has expanded in space and is thus 
sufficiently monochromatic. They grow to an amplitude that in- 
troduces oscillations of the mean velocity of the ions of up to 
c/5. 

The interplay of the short-scale charge density waves and 
magnetowaves causes the electrons to be accelerated to highly 
relativistic speeds upstream of the forming shock. Similar elec- 
tron acceleration (injection) mechanisms upstream of shocks in- 
volving whistler waves have been pro posed by Levinson ( 1992j); 
iKuramitsu & Krasnoselskikhl (l2005alibl) . The strongest electron 
acceleration is, however, observed at the location where the 
shock-reflected ion beam is developing. The electrons are ac- 
celerated to a peak Lorentz factor of 120 and their energy gain 
is thus comparable to the energy associated with the velocity 
change of the shock-reflected ions. If the electron acceleration 
is accomplished by the electromagnetic fields that reflect the in- 
coming upstream ions, then the energy gain of the electrons may 
scale with the ion mass. We may expect in this case that the elec- 
trons are accelerated to a Lorentz factor jm ~ 120OTp/m, that is 

~ 550 if we would use the correct proton to electron mass 
ratio. 

It is interesting to see if this type of electron acceleration 
can also occur close to Solar system shocks. Let us consider 
the Earth bow shock as one of the best known collisionless 
shocks and let us assume a Solar wind speed of 4 x 10^ m/s 
to 7.5 X 10^ m/s. A specular reflection of the incoming Solar 
wind protons by the bow shock would change their energy by 
about 0.8 keV to 3 keV. Electrons with such energies are ob- 
served in a thin sheet close to the shock surface of perpendicular 
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shocks d Anderson et al.|[T979l) . Perpen dicular shocks are cap able 
to produce shock-reflected ion beams ('Lembeg e et al.l2004l) and 
the electron acceleration mechanism we observe here may work 
also at the Earth bow shock. 

4.5. Future work 

This simulation study was concerned with the collision of two 
plasma clouds at a mildly relativistic speed. Its purpose has been 
to better understand the conditions and the mechanisms involved 
in the formation of a shock. This shock will move at an essen- 
tially nonrelativistic speed below the initial collision speed and, 
thus, be relevant for fast SNR flows. A two-dimensional sim- 
ulation geometry was necessary to assess the importance of the 
multi-dimensional filamentation instability for the shock dynam- 
ics. A quasi-parallel guiding magnetic field was used to slow 
down this filamentation, resulting in a planar (one-dimensional) 
shock. 

The formation of the shock could be observed, but the sim- 
ulation had to be stopped at this time due to computational con- 
straints. This simulation has, however, revealed several aspects 
that should be examined in more detail in more speciaUsed sim- 
ulation studies. 

The EES probably formed in response to our initial con- 
ditions, i.e. the sharp jump in the convection electric field at 
the cloud collision boundary. It has to be investigated if the 
EES also forms if this jump is decreased, for example by a 
higher-order field interpolation scheme or by a gradual de- 
crease of the convection electric field by a smo oth change in the 
plasma convection sp eed. The simulations by Bessh o & Ohsawal 
( Il999l) : IZindo et all (12005 ) suggest that this will leave un- 
changed the magnetic field amplification and the electron accel- 
eration. However, the EES may not be so strong and coherent. 

The magnetic field amplitude in the present study is higher 
than it is realistic for a SNR scenario. Future studies must ad- 
dress how the shock formation depends on lower amplitudes of 
the guiding magnetic field. Computationally inexpensive para- 
metric simulation studies that resolve only one spatial dimension 
may provide insight. Initial studies not discussed here indicate 
that the shock formation is delayed by a decreasing magnetic 
field amplitude. 

It is also necessary to follow the plasma colUsion for a longer 
time. An important aspect is here how far the EES can expand 
upstream and how strong the downstream magnetic field is. The 
planarity of the EES and of the shock boundary may permit us 
to use one-dimensional simulations, by which we can expand by 
at least an order of magnitude the box size along the collision 
direction. 

An one-dimensional simulation also allows us to examine 
with a larger number of particles per cell and, thus, lower noise 
levels the secondary instabilities driven by the EES. We have 
found evidence for an instability of the EES to a four- wave inter- 
action. Lower noise levels would allow us to compare the ampli- 
tudes and phases of the EES with those of the secondary waves, 
which is necessary to demonstrate a coherent wave interaction. 
Extending the simulation time together with suitable initial con- 
ditions would, potentially, allow us to investigate what happens 
if the speed of the EES decreases below the local Alfven speed. 
It is possible that the EES decouples from the shock and propa- 
gates independently in form of an Alfven wave packet. 

Finally, it would be interesting to reduce the collision speed 
to about c/10, which is close to the expansion speed of the SNR 
shock, to see if and how many electrons are accelerated to rel- 
ativistic speeds. This will provide insight into the electron in- 



jection efficiency of obUque shocks and, thus, into the abiUty of 
SNR shocks to accelerate electrons to cosmic ray energies. 
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